function q1=q1(x)

q1=2*(3*x(1,1)^2+2*x(1,1)*x(1,2)+x(1,2)^2)/sqrt(x(1,1)^2+x(1,2)^2);